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A new theory is proposed to interpret nuclear spin- lattice relaxation-time (Ti) measurements 
on the spin-1 quasi-one-dimensional Heisenberg antilerromagnet Ni(C2HgN2)2N02(C104) (NENP). 
While Sagi and Affleck pioneeringly discussed this subject in terms of field-theoretical languages, 
there is no theoretical attempt yet to explicitly simulate the novel observations of T^^ reported by 
Fujiwara et ai. By means of modified spin waves, we solve the minimum of T^^ as a function of 
an applied field, pending for the past decade. 



PACS numbers: 76.60.-k, 76.50.-hg, 75.10.Jm 

Predicting a striking contrast between integer- and 
half-odd-integer-spin one-dimensional Heisenberg anti- 
ferromagnets, Haldane [1,2] sparked renewed interest 
in low-dimensional quantum magnetism. The Haldane 
gap, that is, a magnetic excitation gap immediately 
above the ground state, was not only calculated by var- 
ious numerical tools [3-8] but indeed observed in spin-1 
quasi-one-dimensional Heisenberg antiferromagnets such 
as CsNiCla [9] and Ni(C2H8N2)2N02(C104) (NENP) 
[10]. The valence-bond-solid model [11,12] significantly 
contributed toward understanding novel features of the 
Haldane massive phase such as fractional spins induced 
on the boundaries [13,14], a string order hidden in the 
ground state [15,16] and magnon excitations against the 
hidden order [17-20]. The nonlinear-cr-model quantum 
field theory [21,22] skillfully visualized the competition 
between massive and massless phases, while a general- 
ized Lieb-Schultz-Mattis theorem [23,24] gave a criterion 
for the gap formation in a magnetic field. 

Recent progress in the experimental studies also de- 
serves special mention. Not only the single-magnon dis- 
persion relation [25] but also the two-magnon contin- 
uum [3,26] was directly observed by inelastic-neutron- 
scattering measurements on NENP [27] and CsNiCla [28]. 
An applied magnetic field may destroy the Haldane gap 
and bring back magnetism to the system. Such a field- 
induced long-range order was indeed realized in a nickel 
compound Ni(C5Hi4N2)2N3(PF6) [29]. Another family 
of linear-chain nickelates of general formula i?2BaNi05 
[R = rare earth or Y) [30] exhibited a novel scenario of 
one- to three-dimensional crossover [31,32]. When the 
nonmagnetic Y^+ ions are substituted by other magnetic 
rare-earth ions in Y2BaNi05 with a disordered ground 
state, there appears a three-dimensional long-range or- 
der, while the one-dimensional gapped excitations persist 
both above and below the Neel temperature [33-35] . 

Nuclear spin-lattice relaxation time (Ti) has also been 
measured on Haldane-gap antiferromagnets. The field 
[H) dependence of T^^ is of particular interest at both 
high and low temperatures. The high-temperature re- 
laxation rate has been discussed in the context of trans- 
port properties. Takigawa et al. [36] measured the high- 



temperature relaxation rate of ^^P and ^'^V nuclei of 
AgVP2S6, which is also an ideal spin-1 Haldane-gap 
antiferromagnet, and observed the diffusive dynamics 
Tf ^ cx [37]. There is a hot argument [38-41] 

whether the spin transport in quantum spin-gapped an- 
tiferromagnets should be diffusive or ballistic at finite 
temperatures. On the other hand, Gaveau et al. [42] and 
Fujiwara et al. [43,44] measured the low-temperature re- 
laxation rate of ^H nuclei of NENP. As an applied field 
increases, the first excited state moves down and then 
crosses the ground-state energy level. Indeed T^^ reaches 
a peak near the critical field Hc = Aq/ g^B at every tem- 
perature, but it is not a monotonically increasing func- 
tion of H, at low temperatures T < A^/kB in particular. 
Sagi and Affleck [45] formulated the nuclear magnetic 
resonance in Haldane-gap antiferromagnets in terms of 
field-theoretical languages. However, few investigations 
have followed their pioneering argument and there is no 
attempt yet to explicitly fit a theory for the above obser- 
vations. 

In such circumstances, we revisit the low-temperature 
nuclear magnetic relaxation in Haldane-gap antiferro- 
magnets with particular emphasis on its field depen- 
dence. Excluding any phenomenological assumption 
from our argument, we calculate the nuclear spin- lattice 
relaxation rate by means of modified spin waves. A 
new theory claims that as an applied field increases, T^^ 
should initially decrease logarithmically and then increase 
exponentially, well explaining experimental observations. 

We employ the spin-1 one-dimensional Heisenberg 
Hamiltonian 

L L 

jY^Si-Si+^-giiBHY^St . (1) 

1=1 1=1 

In order to illuminate the essential relaxation mechanism 
in spin-gapped antiferromagnets as analytically as possi- 
ble, we do not take any anisotropy into consideration in 
this letter. Magnetic anisotropy quantitatively affects the 
gap amplitude but has no qualitative effect on the whole 
scenario. Scaling temperature and an applied field by the 
Haldane gap, we present a universal theory. Quantitative 
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refinement of the final product in the presence of single- 
ion and orthorhombic anisotropy terms will be considered 
elsewhere. Introducing the Holstein-Primakoff bosonic 
operators 



^2n-l ~ '^26' — alittn CLn , 82^-1 — S — (l|jCln : 

^2+ = bi yj2S-hibn , SL = -S + bibn , 



(2) 



and retaining only bilinear terms of them, we rewrite the 
Hamiltonian as 



n = -2JS^N 

+{2JS + gti^H) 4«n + {2JS - qi^bH) ^ b% 



N 



N 



n=l 



N 

+JSY (albl + Unbn + blal^_^_-^ + 6„a„+i^ , (3) 



n=l 



where N = L/2. In order to preserve the up-down sym- 
metry, or the sublatticc symmetry, we optimize the spin- 
wave distribution functions constraining the total stag- 
gered magnetization to be zero [46-49]: 



Y{aian + bibn)=2NS. 



(4) 



The constraint is enforced by introducing a Lagrange 
multiplier and diagonalizing an effective Hamiltonian 



n = n + 2j\J2 {aUn + bibn) . 



(5) 



Via the Bogoliubov transformation 



_L^e''=(2"-i/2)a^ = afeCosh6»fe - ^sinh^fc , 



_^^g-ife(2n+i/2)^^ = I3kcosh0k-alsmhek, 



(6) 



with tanh26'fe = Scosk/{S + X), we reach the spin- wave 
Hamiltonian 

H = -2JS^N - 2J{S + X)N+jY,^k 

k 

+J E {^k(^W + ^t^k^k) , (7) 

k 

where 



uj^ T gi^nH/J = 2^{S + A)2 - 52cos2A; = Uk • (8) 

Minimization of the free energy gives the optimum dis- 
tribution functions as = (e'^'^fe Z*^^^ — and A is 
then determined through 

Y^i^k +nt + l)cosh20k = {2S+l)N. (9) 



Now we calculate the relaxation rate in terms of the 
modified spin waves. Considering the electronic- nuclear 
energy-conservation requirement, the Raman scattering 
predominates in spin-gapped antiferromagnets. The Ra- 
man relaxation rate is generally given by 



J_ _ 47r(.g^B?i7N)^ 



-Ei/kBT 



hi 



-Ei/kTsT 



(10) 



where A; is the dipolar coupling constants between the 
nuclear and electronic spins in the Ith site, = 7n^ is 
the Larmor frequency of the nuclei with 7n being the gy- 
romagnetic ratio, and the summation ^ . is taken over all 
the electronic eigenstates \i) with energy Ei. By means 
of the modified spin waves, eq. (10) is rewritten as 
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47r(ffMB?t7N)^ 



6{Ej - E,- Hlon) 



0'lEfc,fc'^fe'-fe[("l'"fe - /3j,/3fe)cosh6lfe'Cosh6lfe 

-{ak'al - f3k'pi)smh9k'smh9k] \i) ^ , (11) 



where Aq = Y^i&'^'-Ai. The Fourier components of 
the hyperfine coupling constant exhibit little momen- 
tum dependence when the nuclei take unsymmetrical po- 
sitions in the crystal, which is the case with the pro- 
tons in NENP [44]. Hence we assume in the following 



that Aa 



A 



q=Q 



A. Due to the significant differ- 



ence between the electronic and nuclear energy scales 
(TiwN ^ 10"^ J), oq. (11) ends in 



J_ ^ 4(ff/XB?i7NA)^ 



7r/2 



7r/2 y/v'^k'^ + vhujf^ 



dk, (12) 



where assuming moderate temperatures ksT <C 
^k=-K/2 ~ 2</(S' + A) — gusH, we have approximated the 
dispersion relations as 



J(j^ A + vk^ ± qhbH , 



with 



A = 2J^\{2S + A) , 



V = 



VA(25 + A) ■ 



(13) 



(14) 



Equation (12) claims that an applied field produces two 
distinct effects on Tj~^, one of which originates from the 
Zeeman energy and appears in n^, while the other of 
which appears via the nuclear spins giving the charac- 
teristic weight (u^fc^ + vhcv-^)^^^"^ to the electronic tran- 
sition rate n'^{n'^ + !)• The field effect on the nuclear 
spins escapes observation in critical spin chains with 
a linear dispersion at small momenta. The prefactor 
{v'^k'^ + v?ujjn)~^^'^ is the consequence of the nature of 
the delta function, 
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6{x 



(15) 



where Xi is a zero point of an arbitrary regular func- 
tion ,f{x), and therefore generally arises from quadratic 
dispersion relations of the relevant electronic excitations, 
which are the case with ferromagnets [50-52] as well as 
spin-gapped antiferromagnets. 

Let us fit cq. (12) for the proton spin-lattice relaxation- 
time measurements on NENP [44]. Although A and v, 
given in eq. (14), depend on temperature in principle, 
here we fix {v'^k'^ + vhu!N)~^^^ to its zero-temperature 
value in the integration (12), which is well justified for 
fce^r < Z\ — gjJLB,H and allows us to inquire further into 
eq. (12) analytically. Wc compare the calculations with 
the observations in Fig. 1. Assuming that g = 2 and 
J/ks = 55K [44], we have set A equal to 0.024A~^, 
which suggests the distance between the interacting pro- 
ton and electron spins being about 3.5 A and is consis- 
tent very well with the structural analysis [53]. Under 
the present parametrization, the lowest excitation gap is 
given by A(T = 0)/kB = ^o/^b — 4.0 K, which is some- 
what smaller than that of NENP, 12.8 K [44]. However, 
the scaled function A{T)/Ao well reproduces the upward 
behavior of the Haldane-gap mode as a function of tem- 
perature [49]. 

With increasing field, the relaxation rate first decreases 
moderately and then increases much more rapidly, at 
low temperatures in particular. Although we cannot cal- 
culate beyond the critical field = Aq/qiib — 9.5 T 
on the present formulation, our theory well reproduces 
the observations over a wide field range. For ksT < 
A — gunH, the distribution functions may be approxi- 
mated as n^(n^ + 1) ~ 0''^^/''^'^ and therefore eq. (12) 
can further be calculated as 



1 



8(g/XB^7N^)^ 
irJiv 



-^""■^»K^)-»(S&). 

(16) 



where Kq is the modified Bessel function of the second 
kind. Provided hw^ <C ksT < Z\ — g^sH, we further 
reach 



1 ^ 8((?/XB?l7N>l)'g_Ao/fcBTpQgjjJ^ffMB^f\ 

0.80908 -Inf^) 



(17) 



Thus, as H increases, T^^ should initially decrease log- 
arithmically and then increase exponentially. Figure 2 
suggests that the initial decrease of Tj~^ turns from 
T^^ {X —\nH to T^^ cx H^^^^ with decreasing tem- 
perature. The momentum distribution functions are 
peaked at fc = and their peaks arc sharpened as T 



decreases. 



behaves as S{k) in the low-temperature 



limit. When we replace by 6{k), eq. (12) gives a 
if~^/^-linear field dependence of T{~^. 



The temperature dependence is mainly described by 
the term e-{^o-giJ.BH)/kBT ^fut further decorated due 
to the temperature-dependent energy spectrum. The 
inelastic-neutron-scattering peak position of the lowest- 
energy mode exhibits an upward behavior with increas- 
ing temperature, for fceT > Ao/2 in particular [54 56], 
where the slope of InTj^^ to correspondingly in- 

creases with increasing temperature. 

The nuclear magnetic relaxation in the Haldane-gap 

antiferromagnet NENP has been interpreted in terms of 
a modified spin-wave theory. The field dependence of 
T]~^ was analyzed in detail and the minimum ofT^^ as 
a function of H, pending for the past decade, was solved. 
We consider that such a field dependence of Tj~^ is qual- 
itatively common to spin-gapped antiferromagnets. We 
encourage low-temperature Ti measurements on related 
materials such as the ferromagnetic-antiferromagnetic 
bond-alternating compoimd (CH3)2CHNH3CuCl3 [57] 
and the two-leg ladder antiferromagnet SrCu203 [58]. 
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FIG. 1. Dependences of the nuclear spin-lattice relaxation 
rate on a field parallel to the chain (a) and temperature 
(b) in NENP (symbols) [44] are compared with the modified 
spin- wave calculations (lines), where the field and tempera- 
ture are scaled by the lowest excitation gap Z\o- 
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FIG. 2. The momentum distribution functions (dotted 
lines) and (solid lines). 
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